Charge accumulation on a Luttinger liquid 
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The average charge Q on a quantum wire, modeled as a single-channel Luttinger liquid (LL), connected to 
metallic leads and coupled to a gate is studied theoretically. We find that the behavior of the charge as the 
gate voltage Vg varies depends strongly on experimentally adjustable parameters (length, contact transmission, 
temperature,. . . ). When the intrinsic backscattering at the contacts is weak (i.e. the conductance is close to 
jh at high temperature), we predict that this behavior should be described by a universal function. For short 
such wires, the charge increases roughly linearly with Vg, with small oscillations due to quantum interference 
between electrons scattered at the contacts. For longer wires at low temperature. Coulomb blockade behavior 
sets in, and the charge increases in steps. In both limits dQ / dVc which should characterize the linear response 
conductance, exhibits periodic peaks in Va- We show that due to Coulomb interactions the period in the former 
limit is twice that of the latter, and describe the evolution of the peaks through this crossover. The study can 
be generalized to multi-channel LLs, and may explain qualitatively the recent observation by Liang et aim of a 
four-electron periodicity for electron addition in single-walled carbon nanotubes. 

PACS numbers: 71.10.Pm,73.23.-b,73.23.Hk 



I. INTRODUCTION 

The conductance of metallic single-walled carbon nan- 
otubes has been shown to depend strongly on the nature of 
the contacts between the nanotube and the leads. In a typi- 
cal experimental setup a bias voltage is applied across a nan- 
otube connected to metallic leads, while a gate voltage applied 
to a third electrode acts as a chemical potential and modu- 
lates the charge on the nanotube. ' --^ ^.s wjjgjj ^jje contacts be- 
tween the nanotube and the leads are poor. Coulomb block- 
ade behavior sets in, and the conductance exhibits a series of 
sharp peaks as the gate voltage increases. ^-'•^ In contrast, the 
conductance of devices with near-perfect contacts is close to 
the theoretical maximum of Ae'^/h for all gate voltages, with 
small quasi-periodic oscillations due to Fabry-Perot electron 
interference.^ 

An interesting question that arises from these experiments 
is how the conductance of the system evolves in between 
these limits. Simple arguments indicate that the conduc- 
tance undergoes a change in periodicity between the Coulomb 
blockade regime and the Fabry-Perot limit. Deep in the 
Fabry-Perot limit, based on a non-interacting picture in which 
each conducting channel is approximately independent, one 
expects one electron per channel (i.e. 4 for a nanotube) 
to be added to the wire per period of conductance oscil- 
lation. In the Coulomb blockade limit, peaks occur upon 
each electron addition process. This reasoning is in accord 
with recent experimentsif^, in which transport measurements 
were performed on devices exhibiting a broad range of low 
temperature conductances. These experiments confirm the 
Fabry-Perot picture in the most conducting samples, regular 
Coulomb blockade behavior in the least conducting ones, and 
observed an interesting clustering of peaks into groups of four 
in an intermediate limit'. 

Theoretically, there are two qualitative issues brought up 



by these experiments. First, why should the simple Fabry- 
Perot type structure obtained from a non-interacting quasipar- 
ticle scattering approach apply even in the highly conducting 
samples? This experimental result is somewhat surprising, 
since both theoretical expectation and numerous experiments 
indicate that Coulomb effects in nanotubes modify them from 
Fermi to Luttinger Liquids, and in particular should not have 
electron-like quasiparticles! This question was addressed al- 
ready in Ref. 7, where it was demonstrated that, for highly 
conducting samples, the behavior expected from LL theory (at 
least at low source-drain bias) is qualitatively indistinguish- 
able from the naive quasiparticle scattering prediction. With 
this understood, the second issue raised is the nature of the 
crossover from this Fabry-Perot limit to the Coulomb block- 
ade. Qualitatively, one would like to understand how the po- 
sitions of the peaks evolve between the two regimes. This key 
qualitative issue of how the center of the peaks evolve is the 
focus of our work. Of course, the full universal crossover is 
interesting, but is more quantitative than qualitative. 

A proper theoretical treatment of this evolution would re- 
quire calculating the conductance of nanotube devices with 
arbitrary contact resistances. While this has been done pertur- 
batively in the case of near-perfect contacts^, non-perturbative 
techniques are required for intermediate contact resistances. 
Finding the conductance in the crossover regime is therefore 
of considerable difficulty. For low bias voltages, however, the 
qualitative features of the conductance should be manifested 
in the derivative of the average charge on the nanotube with 
respect to the gate voltage. This can be determined from an 
equilibrium calculation of the free energy and is consequently 
more tractable. Even with this simplification the problem is 
nontrivial, so we use Feynman's variational principle**-' to cal- 
culate the free energy and subsequently the charge. Since 
there is both experimental"* and theoreticali2iii evidence that 
carbon nanotubes behave as Luttinger liquids '^(LLs), we con- 



sider more generally the average charge on a single-channel 
LL (with spin) with arbitrary electron-electron interaction 
strength. Our results can be generalized to multi-channel LLs 
and therefore can be applied to nanotubes, which have two 
channels of conduction near the Fermi energy. 

As expectedi^ in the Coulomb blockade regime, when the 
contact resistance is large, we find that the charge increases 
discontinuously in steps whenever the gate voltage (in units in 
which two electrons are added to the wire whenever Vq is in- 
creased by one unit) Vg = (2n + l)/4, where n is an integer 
The discontinuous character of these charge jumps is an arti- 
fact of our variational method and translates into (5-function 
peaks in cr = dQ/dVa- The charge should be a continuous 
function of Vg, and the height of the peaks in a should be fi- 
nite. The variational method predicts accurately the locations 
of these peaks, but not their heights. If we correct for this ar- 
tifact, cr should behave qualitatively as depicted in curve 1 in 
Figure^ Despite the spin degeneracy, each peak corresponds 
to the addition of a single electron to the wire. Indeed, in this 
case single electrons can enter the LL only when they acquire 
the additional energy required to overcome the Coulomb re- 
pulsion energy and in this regime the periodicity of the peaks 
is 5Vg = 1/2. 

Although the method we use has the advantage of allowing 
us to consider arbitrary contact resistances, it yields unphysi- 
cal results when the contact resistance is low. In this regime 
we resort instead to perturbation theory . In this limit we find 
that the effects of the electron-electron interactions are small. 
The spin up and down electrons propagate independently of 
each other so that two electrons are added to the LL per Fabry- 
Perot oscillation in a. Consequently, a exhibits a broad sinu- 
soidal oscillation'** with a period Vg = 1, which is twice as 
large as the period in the Coulomb blockade regime. This be- 
havior is represented by the curve 4 in Fig. 1 . 

The aim of this paper is to address how the crossover be- 
tween the Coulomb blockade and the nearly-perfect contacts 
regime takes place. One possibility is that as the contact re- 
sistance decreases, the Coulomb blockade peaks remain fixed 
in position but broaden asymmetrically so that they eventually 
combine to form the Fabry-Perot oscillations. A second pos- 
sibility is that the Coulomb blockade peaks shift in position 
as they broaden and at some point collapse into each other to 
form the Fabry-Perot peaks. We find the latter possibility to be 
the case, and we illustrate it schematically in Figure^ which 
contains a sketch of cr = dQ / dVG (or the conductance G) as a 
function of Vg for different contact resistances. As the contact 
resistance decreases, we find that the peaks begin to shift to- 
ward the nearest half-integer value of Vq and the magnitude of 
the charge jump diminishes. Correcting for the discontinuity 
in the charge, a should evolve qualitatively as shown in Figure 
n (curves 2 and 3) when the contact resistance decreases. 

It should be stressed that in the limit of weak backscat- 
tering at the contacts, this crossover is universal. This fol- 
lows from general renormalization group reasoning. It is well 
known that, for a spin- 1/2 LL with full spin-rotational in- 
variance (and this generalizes to the LL-Fermi liquid contacts 
considered here), there is a single relevant impurity backscat- 
tering operator. Hence, if the strength of the backscattering 
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a = dQ/d% 
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FIG. 1: The qualitative features of the crossover between good and 
poor contacts. The derivative of the charge with respect to the gate 
voltage (or the conductance) is depicted schematically for various 
values of the backscattering strength at the contacts. Curve 1 rep- 
resents the Coulomb blockade situation, and curves 2 and 3 depict 
progressively smaller backscattering. Curve 4 illustrates the Fabry- 
Perot oscillations present in the low backscattering limit. 

at each contact is weak and they are well-separated from one 
another, all other effects of the imperfections of the contacts 
scale rapidly to zero on the scale of wire length. Hence, the 
behavior of such samples should be well-described by a model 
containing only the ideal LL Hamiltonian and this leading 
backscattering operator at each contact. In particular, then, 
the charge Q on the wire is a universal function of the in- 
teraction strength g, the renormalized backscattering strength 
^ uL^^^, where L is the length of the wire, the gate voltage 
Vg, the thermal length kBTL/hvp, etc. This still contains 
many parameters that will affect the quantitative form of this 
crossover, but we expect our results for the evolution of the 
peak positions to be much more robust. 

We establish the validity of our technique in Section II, 
where we apply Feynman's variational principle to calculate 
the free energy of an infinite, spinless LL with a single impu- 
rity. In Section III A, we discuss the variational solution for 
the charge on a spinless LL connected to semi-infinite leads. 
A similar analysis is performed in Section III B where we take 
into account spin. Finally, our results are summarized in Sec- 
tion IV. 



11. INFINITE LUTTINGER LIQUID WITH A SINGLE 
IMPURITY 

In this section we use Feynman's variational principle to 
calculate the free energy for an infinite, spinless Luttinger 
with an interaction parameter g describing the strength of 
the interactions and a single impurity of scattering strength 
u. Although the free energy for this system can be computed 
exactly''', applying the variational principle here allows us to 
establish the validity of this technique by comparing our re- 
sults to the well-known RG scalings resultsi^*ii. 

We begin by writing an effective bosonized actioni^ for this 
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system, 



5 = /jV— |eoK)|'-u r drcos[2eoiT)], (1) 
^ ^9 Jo 



where represents 9{x) evaluated at the impurity site. The 
bosonic field 6 and its dual (j> are related to the right- and left- 
moving fermionic fields 4'r/l via the transformation ipn/L ~ 
gi(0±e) 'pjjg Fourier transform conventions we use through- 
out are 



1 



= 



(2) 



We use non-standard conventions for simplicity so that both 
0{uj) and 6{t) are dimensionless. 

Determining the free energy for this system is nontrivial 
due to the presence of the cosine term in the effective action. 
We therefore approximate the free energy using Feynman's 
variational principle^'^, which states that for some trial action 
S", the exact free energy obeys the following inequality, 



F < 



1 



InZ' + -{S-S')i 



F„ 



(3) 



where Z' is the partition function corresponding to S" {Z' = 
J VOe^^ ). The trial action S' is chosen to have a tractable 
form and to depend on variational parameters that are deter- 
mined by minimizing the variational free energy, F„. The re- 
sulting variational free energy will be the best estimate for the 
exact free energy given the form of the trial action used. 

In our case, the most general tractable trial action is 
quadratic in 6, so we take our trial action to be 



Tig 



(4) 



where A is a non-negative, frequency-independent variational 
parameter This form can be obtained by expanding the cosine 
term in the effective action to second order in 9 and replacing 
the scattering strength u by an effective scattering strength 
A/2. This turns out to be the most general quadratic action 
that one needs to consider Even if one takes into account an 
explicit frequency dependence for A, the solution that mini- 
mizes the free energy is frequency-independent. 

Computing the variational free energy using Eq. @ and set- 
ting dFy/d\ — 0, we find in the zero-temperature limit that 
the free energy is minimized when A satisfies the following 
equation. 



\ = 2u 



A 



A + eo/(7rg) 



(5) 



where eg is a high-energy cutoff. For g > 1, Eq. Q is satisfied 
only when A = 0. Thus the variational free energy is min- 
imized by choosing a trial action which completely neglects 
the effects of scattering. For g > I this result is consistent 



with renormalization group arguments'*, which state that u is 
irrelevant for the case of attractive interactions and that the 
scattering strength is therefore renormalized to zero at ener- 
gies much smaller than eg- However, when g — 1 the system 
reduces to a one-dimensional free fermion gas impinging on 
a barrier of height u. In this case the scattering strength is 
marginal, and the physics of the system (the value of the con- 
ductance, etc.) should depend on the the scattering strength, 
which our variational technique fails to predict. For g < 1, 
there exists a nontrivial solution to Eq. (|5) given by 



A « 2u 



27r.guN^ 9/(1-9) 



(6) 



that minimizes the free energy. It follows from Eq. ^ that A 
increases from at g = 1 to 2u at g = so that the effective 
scattering strength increases with the strength of the interac- 
tions. This behavior is also consistent with renormalization 
group arguments since u is relevant in this range of g. For this 
simple system, our variational technique therefore reproduces 
the well-known scaling results'* except in the case 5 = 1. 



III. FINITE-SIZE LUTTINGER LIQUID CONNECTED TO 
LEADS AND COUPLED TO A GATE 



A. Spinless fermions 

In this section we calculate the charge on a spinless LL of 
length 2L and interaction parameter g connected to two semi- 
infinite leads characterized by an interaction parameter g^. 
For simplicity, we assume that the Fermi velocity v is uni- 
form in the LL and the leads. The contacts with the leads at 
positions x — ±L are modeled as impurities of equal scat- 
tering strength u, measured in units of eL = v/2nLg^. A 
gate voltage Vq is applied to the LL, where Vq is expressed 
in units of tt^el so that one electron is added to the LL when 

Vg^Vg + I. 

The effective action for this system is 



S 



EE 

uj-n, a— lb 

^ f drY, cos 

-^0 C=±i 



6a{uJn)\'^ICa{iOn ) - v^nVcO-iO) 

V2e+{T)+CV29-{T))],i7) 



where e± ^ [e{x = L) ± 9{x = -L)]/V2 and 9^{0) repre- 
sents 9-{uJn = 0). The functions /C±(cj) appearing in Eq. 
are given by 



27rer 



- 1 
■9L 



1 r 

9 



tanh ( - 



-1 



V27reL5/ 



(8) 



To calculate the variational free energy, we take a trial action 
of the form 



V27r[(VG + AG)f?-(O) + /i+0+( 



(9) 
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and assume that the variational parameters A± are frequency- 
independent and non-negative. If one allows A± to depend on 
frequency, one finds as in the single-impurity case that the free 
energy is minimized when X± are frequency-independent. 
The quadratic terms in Eq. (|9j can be obtained by expand- 
ing the cosine term in the effective action to second order and 
replacing u by A±/2, so we can once more interpret X± as ef- 
fective scattering strengths. To obtain the linear terms in Eq. 
(|9}, we shift Vq by a variational parameter and introduce 
another variational parameter /i+ that multiplies 6*^(0). The 
inclusion of the latter term, which does not appear in the ef- 
fective action in Eq. Q, is necessary in order to preserve the 
invariance of the original Hamiltonian under the addition of 
an extra unit of charge. 

The variational free energy computed from the trial action 
in Eq. (|9} is given in Eq. iAl\ in Appendix A. Setting the 
derivatives of the free energy with respect to the variational 
parameters to zero, we find that the free energy will be mini- 
mized when A± satisfy the following equation. 



A+ = A_ = A = CtA'^g 



Ap-/(A) 



X cos ttVg + jyVO^A^^^F^^W^A^ 



(10) 



where cos(7r/x_|./A) = ( — ±1, 77 — ±1, and A = 
"^dLd/ {9l + <?)■ The dimensionless parameter that determines 
the behavior of the system is 7 = 2u{nAeL/eQ)^, and in the 
zero-temperature limit /(A) is given by 



/(A) 



dx 



a=± 



K.ai2TTeLx) + X x/A + X/2 



.(11) 



The signs of ( and 77 depend on the gate voltage and should be 
chosen such that the free energy is minimized and A is non- 
negative. To satisfy these conditions, we choose ( ^ +1 if 
— 1/2 < Vq < 1/2 and change the sign of C whenever Vq 
Vg + 1. Also, we take ry = — 1 if < Vg < 1/2 and change 
the sign of 77 whenever Vq — > Vq + 1/2. 

Since the gate voltage acts as a chemical potential, the av- 
erage charge on the LL (in units of the electron charge) is 
Q = — (9i^-u/(9VG)/(7r^ei). When the free energy is mini- 
mized, the charge can be written in terms of A as 



QW 



Vg + Xc 
1 + A 



VG + -v/7^A2Ae-2/(A) _A2. (12) 

TT 



To determine the charge, one must therefore find the solution 
to Eq. dlOt that minimizes the free energy. 

In the limit 7 1 we can solve Eq. ilO\ analytically 
by assuming that A is large and roughly independent of Vq- 
The dominant contribution to the integral in /(A) then comes 
from the region where x ^ 1. We therefore approximate 
K.±{2'neLx) « 2a;/A and take /(A) w 0. With these assump- 
tions we find that A = is the only solution when A > 1. 
For A < 1, which is the relevant physical situation since we 
are interested in repulsive interactions, the free energy is min- 
imized by a nontrivial value of A. The charge is then 



n + rj 



1 



1/(1-A) 



(13) 



where n is the closest integer to Vg ■ This solution corresponds 
to the Coulomb blockade limit since the charge increases in 
steps whenever Vg equals a half-integer. We note that the 
charge jumps discontinuously at these values of Vq, resulting 
in the emergence of J-function peaks in a. This discontinuity 
is an artifact of the variational method and can be traced back 
to the replacement of the cosine term in the effective action 
by quadratic terms in our trial action. Intuitively, the charge 
should be a continuous function of Vg, and the peaks in a 
should be rounded and have finite height. 
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FIG. 2: The charge Q and a — dQ/dVc as functions of the gate 
voltage for a system with gL ~ 1, g ~ 0.25, el/eo ~ 2.55 x 10~^, 
and 7 — 2.37. The discontinuity in the charge at Vb = 1/2 is 
an artifact of the variational technique that arises from our replacing 
the cosine term in the effective action by quadratic terms in our trial 
action. This results in 5-function peaks in a which should instead be 
rounded and have finite height. 
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FIG. 3: The charge Q and a — dQ/dVa as functions of the gate 
voltage for a system with gL ~ 1, g ~ 0.25, ei /eo = 2.55 x 10~^, 
and 7 — 0.8. As 7 decreases the discontinuous jump in the charge 
at Vg = 1/2 diminishes, and the slope of the charge away from 
Vg = 1/2 increases. 

For smaller values of 7, we obtain the charge by numeri- 
cally solving Eq. ilQ\ for the value of A that minimizes the 
free energy. To illustrate the evolution of the charge, we 
plot Q and its derivative a = OQ/OVg as functions of Vg 
for different values of 7, focusing on a system at zero tem- 
perature with Fermi liquid leads {g^ ~ 1), g = 0.25, and 
fi/eo = 2.55 X 10"'^. The solid lines in Figures|2]and|3lrep- 
resent Q and a evaluated at 7 = 2.37 and 0.8, respectively. 
For comparison, these quantities evaluated at 7 = are shown 
as the dashed lines. These figures illustrate that the magnitude 
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of the discontinuous jump in the charge, 5Q, diminishes as 7 
decreases. Additionally, the slope of the charge away from 
half-integer values of Vq increases from toward 1 . 

As 7 decreases further, the behavior of the variational solu- 
tion for the charge depends on whether gj^ is greater than or 
less than 1/2. For < 1/2, 6Q decreases smoothly to zero 
as 7 — > 0. This is illustrated in Figure |4lb), which contains 
6Q versus 7 when g^ — 1/4. Apart from the discontinuity in 
the charge at half-integer values of Vg, in the low backscatter- 
ing limit the variational method produces the expected Fabry- 
Perot oscillations. 
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FIG. 4: The magnitude 5Q of the discontinuous jump in the charge 
at half-integer values of Vg as a function of 7 for a system with 
e^/eo = 2.55 x 10"^ and a) ql ^ I and b) gL = 0.25. 

With g^ > 1/2, SQ drops abruptly to zero when 7 equals a 
critical value denoted by jc- This feature is illustrated in Fig- 
ure |4l a), where 7c « 0.78. Below 7c, the variational method 
predicts unphysical behavior in the charge that, in particular, 
is inconsistent with a perturbative calculation of the charge in 
the limit 7 <C 1. Details of the variational solution in this 
regime will therefore be deferred to Appendix B. The varia- 
tional method predicts the presence of two sharp peaks in a 
per period rather than the single broad peak expected from 
the Fabry-Perot oscillations. Figure |5] contains Q and cr for 
a system with g^ ~ 1 and 7 < 7c. The double-peak struc- 
ture that emerges is a consequence of the free energy being 
minimized when A = in the region between these peaks. 
We attribute this shortcoming of the variational technique to 
the method failing to capture the analytic terms in the free en- 
ergy, which presumably dominate in this regime. In contrast, 
when gL < 1/2 even the lowest-order term in perturbation 
theory diverges, so there are no analytic terms in the free en- 
ergy. This explains why the variational technique yields rea- 
sonable results when g^ < 1/2 for arbitrary 7 but fails in the 
low backscattering limit when g^ > 1/2. 

In the limit 7 <C 1, the charge can instead be calculated 
perturbatively in u when 5^ > 1/2. To second order in u, we 
find 

Q^VG-Aj^sm{2TTVG), (14) 

where A is a positive constant. It follows that a is roughly 
constant, exhibiting small Fabry-Perot oscillations with peri- 
odicity SVg = 1. Thus, although the variational method fails 
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FIG. 5: The charge Q and a = dQ/dVc as functions of the gate 
voltage for a system with Fermi liquid leads {ql = 1), g ~ 0.25, 
el/eo = 2.55 X 10""^, and 7 — 0.47 < 7c. The two peaks present 
in a arise from the free energy being minimized when A = in the 
region between the peaks. This double-peak structure is an artifact 
due to the variational method failing to capture the analytic terms in 
the free energy. 



in the low backscattering limit, the expected behavior is cap- 
tured perturbatively. 

In summary, we have shown that for a spinless, single- 
channel LL the periodicity of a is the same in the low 
backscattering and Coulomb blockade limits. After correcting 
for the artifacts of the variational technique, we find that the 
Coulomb blockade peaks in a simply broaden symmetrically 
as the backscattering strength decreases, eventually evolving 
into the broad Fabry-Perot peaks in the low backscattering 
limit. 



B. Spinful fermions 

We now repeat the calculations of Section III A taking 
into account spin. It is convenient to decompose the system 
into charge and spin sectors via a change of basis to fields 
^p/(T = (^T ^ ^i)/v^- Since the interactions only affect the 
charge sector, the spin sector in both the leads and in the LL 
is characterized by an interaction parameter of g'^ = 1. The 
charge sector is characterized by an interaction parameter g^ 
in the leads and g in the LL. 

When spin is taken into account, the effective action in Eq. 
is generalized to 

^ = V |0,K)|2/C,K)-2^VbV(O) 

- ^ Trfr ^cos[0p+(T)+CiVM 

+ C20.+ (r)+CiC2e.-(T)], (15) 

where a = p±, cr± are charge and spin indices respectively, 
and Ci_2 = ±1. The functions ICa{oJn) are defined by 
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where = sl, g"^ = g, gl'^ = g"^ = 1, Cp/.4 
Cp/cr- = —1- We assume a trial action of the form 



1, and 



(0) 



(0)+/i.-^.-(0)].(17) 

Here, \p/a± represent effective scattering strengths and Xq, 
fj.p+, and fj.a-± are additional variational parameters. The vari- 
ational free energy computed from this trial action is given in 
Eq. iA2\ in Appendix A. As in the spinless case, the free en- 
ergy is minimized when the effective scattering strengths are 
equal, so we define Ap/^i = A. 

As outlined in Appendix A, setting the derivatives of the 
free energy with respect to the variational parameters to zero 
leads to two sets of equations for the variational parameters 
that result in physical solutions for the charge. In the first set, 
the charge is given in terms of A by 



2,7 



Q(A) = 2Vg + ^^y2^A+lg-2/'(A) 
TT 

where A is a solution to 

A = cyA(^+i)/2e-^'(^) cos 



7tQ{X) 



(18) 



(19) 



Here, 7' = 2{tA"i/2(7rAeL/eo)('^+i)/2 is the dimensionless 
parameter that determines the behavior of the system, and 



1 



(27reLw) + A 
1 



w/A + A/2 u + A/2 



(20) 



in the zero temperature limit. The values of 77 and ( are the 
same as in the spinless case. In the second set of equations, 
the charge is determined from the following equation. 



^IQ/2 - Vbl = yA(^+i)/2e-^'(^) I cos (^) 

y2AA+lg-2/'(A) sin2 (^^^ - A2 



where A is a function of the charge via 

A(0) = \n{Q/2 - Vg) tan(^g/2)|. 



(21) 



(22) 



To determine which set of equations leads to an absolute min- 
imum of the free energy, we find numerically the minimum 
free energy solution to both sets and retain the one with the 
lower free energy. 

In the Coulomb blockade regime where 7' 1, the charge 
is given by Eq. ( I18> when n — 1/4 < Vg < n + 1/4 and 
by Eq. J21> when n + 1/4 < Vg < n + 3/4, where n is an 
integer. In this limit the charge can be retrieved analytically 
by assuming A ^ 1 and taking /'(A) w 0. As in the spinless 
case, we find nontrivial solutions for A only when A < 1. The 
charge in this range of A is given by 



n + 2r/ 



n/2 - Vg 



/2/(l-A) 



(23) 



where n is the closest integer to 2Vg and 77' = —1 if < 
Vg < 1/4 and changes sign whenever Vg + 1/4. 

Thus the charge increases in steps, jumping discontinuously 
by an amount 5Q « 1 at half-integer values of 2Vg- This 
discontinuity is again an artifact of our technique and results 
in (5-function peaks in a rather than rounded peaks of finite 
height. These (5-function peaks occur at half-integer values of 
2Vg, so the period of a in the Coulomb blockade regime is 
5Vg = 1/2. 

As 7' decreases and the system moves away from the 
Coulomb blockade limit, the window of Vg inside of which 
the second set of equations yield the lowest free energy begins 
to shrink. The charge remains discontinuous at the endpoints 
of the window, causing the peaks in a to shift. Figures |6j0 
and |8] represent 7' — 1.56, 0.67, and 0.58, respectively, and 
illustrate how Q and a evolve as 7' decreases. For our nu- 
merical data, we again focus on a system at zero temperature 
with Fermi liquid leads, g = 0.25, and /co = 2.55 x 10^"^. 
To illustrate the nonlinear behavior in the charge more clearly, 
Q — 2Vg has been scaled by a factor of 3 in Figure^la) and 
a factor of 4 in FigurelSJa). In Figure|SJ the peaks in a have 
shifted significantly and appear at roughly Vg = 0.4 and 0.6. 
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FIG. 6: The charge Q and a — dQ/dVc as functions of the gate 
voltage for a system with ql ~ 1, g ~ 0.25, eL/eo ~ 2.55 x 10~^, 
and 7' = 1.56. At this value of 7', the peaks in a appear only slightly 
shifted away from Vg = 1/4 and 3/4. 
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FIG. 7: The charge Q and a — dQ/dVa as functions of the gate 
voltage for a system with ql ~ 1, g ~ 0.25, ei/eo = 2.55 x 10~^, 
and 7' = 0.67. For clarity, we have scaled Q — 2Vg by a factor 
of 3 in the plot of the charge. The peaks in a are now located near 
Vg = 0.3 and 0.7. 

When 7' reaches a critical value 7^, the first set of equa- 
tions yields the lowest free energy for all Vg, and the charge 
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FIG. 8: The charge Q and a — dQ/dVc as functions of the gate 
voltage for a system with ql = 1, g ~ 0.25, ^l/^o = 2.55 x 10~^, 
and 7' = 0.58, which is close to 7^. For clarity, Q — 21/0 has 
been scaled by a factor of 4 in the plot of the charge. The peaks in 
o continue to move toward Vg = 1/2 and now appear at roughly 
Vg = 0.4 and 0.6. 
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a) b) 

FIG. 9: The charge Q and a = OQ/OVg as functions of the gate 
voltage for a spinful system with ql = 1, g = 0.25, ei/eo = 2.55 x 
10~^, and 7' = 0.33, which is slightly below 7^. To ejnphasize the 
nonlinear behavior in the charge, we have scaled Q — 2 Vg by a factor 
of 30 in the plot of the charge. The double-peak structure in a is an 
unphysical artifact of the variational method that results from the 
technique neglecting analytic terms in the free energy. 



is given everywhere by Eq. J18> . The discontinuous jump in 
the charge drops smoothly to zero at 7^, and the charge re- 
mains continuous at smaller values of 7'. For 7' < 7^, the 
variational method predicts unphysical behavior in the charge 
that is inconsistent with perturbation theory when 7' <C 1, 
so we defer details of the solution to Appendix B. In contrast 
to the spinless case, where the method was reliable for arbi- 
trary 7' when gL < 1/2, the variational technique fails in this 
range of 7' for all g^. This is consistent with the fact that 
in the spinful case the effective interaction parameter for the 
leads is (g^ + l)/2, which is always greater than 1/2. Figure 
|9] illustrates the unphysical behavior in the charge predicted 
by the variational method when 7' — 0.33 < 7^. For clarity, 
Q — 2Va has been scaled by a factor of 30 in the charge. Non- 
linear corrections to the charge are therefore quite small at this 
value of 7', and a is nearly constant. The variational method 
predicts two peaks in a per period rather than the single broad 
Fabry-Perot peak expected in this regime. This double peak 
structure persists as 7' decreases further until Q — 2Vg and 
the peaks disappear altogether As in the spinless case, we at- 
tribute this failure of the variational technique to the method 
neglecting analytic terms in the free energy. We therefore em- 
phasize instead a perturbative calculation of the charge in the 
low backscattering limit. To second order in u, we find per- 
turbatively that for any g^ > 0, 

Q = 2Vg- sin(2^VG), (24) 

where A' is a positive constant. The period of the Fabry-Perot 
oscillations is therefore 6Vg — 1, which is the same as in the 
spinless case (see Eq. il4\ ). 

In summary, we have shown that the period of a for a spin- 
ful LL increases by a factor of 2 as the system goes from 
the Coulomb blockade regime to the low backscattering limit. 
The numerical data in Figures|6]through|8]along with the per- 
turbative calculation of the charge in Eq. (I24> demonstrate that 
once we coiTect for variational artifacts, this crossover in pe- 
riodicity takes place as sketched schematically in Figure [2 
As the backscattering strength at the contacts decreases, the 
Coulomb blockade peaks shift in position as they broaden. 



and eventually combine to form the broad Fabry-Perot oscilla- 
tions. As mentioned before, the physics behind this crossover 
in periodicity is that in the Coulomb blockade regime elec- 
trons are added one-by-one due to electron-electron interac- 
tions. In the low backscattering limit, however, the spin up 
and spin down electrons essentially propagate independently 
of one another, so that two electrons are added to the LL per 
period of the Fabry-Perot oscillations in a. 



IV. DISCUSSION 

We discussed the average charge on finite-size spinless and 
spinful Luttinger liquids as a function of an applied gate volt- 
age. As expected, in both cases the charge increases almost 
linearly with gate voltage when the contacts are good, while 
the charge increases in steps in the Coulomb blockade regime. 
In the spinless case, the derivative of the charge with respect 
to the gate voltage, cr, has the same periodicity in both lim- 
its. When spin is taken into account, however, we showed 
that each Fabry-Perot peak begins to separate into two peaks 
as the contact resistance increases. As illustrated schemati- 
cally in Figure n the spacing between these peaks increases 
with the contact resistance until the Coulomb blockade limit 
is reached and the period of a is reduced by a factor of 2. The 
physics behind this crossover is that when the contacts are 
good, electrons from the two spin channels essentially prop- 
agate independently of one another so that two electrons are 
added to the LL per period of the Fabry-Perot oscillations. In 
the Coulomb blockade regime, however, despite the spin de- 
generacy electrons are added to the LL one at a time due to 
electron-electron interactions, resulting in a reduction in peri- 
odicity. 

In an obvious generalization, we expect that if there are 
n conduction channels, then the period should change by a 
factor of 2n, with the factor of two arising from spin. This 
crossover in periodicity should also appear in the behavior of 
the conductance. Specifically, rather than pairs of Coulomb 
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blockade peaks combining to form the Fabry-Perot oscilla- 
tions, we expect 2n Coulomb blockade peaks to collapse onto 
each other in the low backscattering limit. Addressing how 
these peaks collapse, however, requires further investigation. 

This crossover in periodicity can be used to interpret re- 
cent experimental data on carbon nanotubes, which have 
been shown to behave as two-channel LLsiiStil. As de- 
scribed in the Introduction, the conductance of nanotubes has 
been measured in devices with near-perfect contacts^, in the 
Coulomb blockade limit^-^, and more recently in an interme- 
diate regime ' . When the contacts are good, we have seen that 
four electrons should be added to the NT per period of the 
Fabry-Perot oscillations in the conductance. As the contact 
resistance increases, corresponding to an increase in 7', each 
of these oscillations should begin to separate into four distinct 
peaks as the crossover begins to take place. The manner in 
which this takes place cannot be implied from our single chan- 
nel calculation, though a few possibilities can occur: all four 
peaks will appear simultaneously and then evolve and separate 
symmetrically; two peaks will arise, evolve and separate sym- 
metrically, then each of them will split into two peaks; two 
peaks will arise and separate first, then when they reached a 
certain distance two more peaks will appear. The recent ob- 
servation of a four-electron periodicity for electron addition 



by Liang et al. may be a manifestation of this crossover (see 
Figure 1 in that reference). The conductance of their devices 
exhibited peaks grouped into clusters of four as the gate volt- 
age increased. We note that their measurements are consis- 
tent with our prediction that the peaks shift in position as the 
contact resistance changes since adjacent conductance peaks 
within a given cluster are closer together than adjacent peaks 
in neighboring clusters. As the contact resistance increases 
further, the four peaks within each cluster should continue to 
separate and become sharper until the Coulomb blockade limit 
is reached and the period of the conductance is reduced by a 
factor of four. 
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APPENDIX A: VARIATIONAL FREE ENERGY 

The variational free energy for the spinless LL computed 
using the trial action in Eq. (|9} is 
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The variational free energy for the spinful LL corresponding 
to the trial action in Eq. ilH is 
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Ap/cr± = A and either 



with a = p±, (T±. Equation iA2\ is minimized when 



sm f ' 



In this case, one can not obtain a decoupled equation for A. 
Instead, we first use the analogue of Eq. iA7\ to write A as a 
function of the charge. We then use the remaining two equa- 
tions to obtain a decoupled equation for the charge, Eq. i21\ . 



= cos (^) = 0. (A4) 



Equation iA3\ leads to the following coupled equations for A, 
/icr_, and the charge Q: 
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By inspection and confirmed by numerical analysis, these 
equations are satisfied when /i^- and A is a solution to 



A = Cl'X^^+'^^'e-''^^^ 
X cos 



ttVg + iWl'^X^+^e-^i'W - A2 



(A8) 



This decoupled equation for A can be used in conjunction with 
Eq. iA5\ to obtain Eq. ( I18> for Q in terms of A. 

Equation JA4> yields a set of coupled equations identical to 
Eqs. ( IA5> -( IX7t except that the cosines are replaced by sines. 



APPENDIX B: VARIATIONAL RESULTS IN THE LOW 
BACKSCATTERING LIMIT 

For completeness, in this Appendix we review the varia- 
tional solution in the limit of near-perfect contacts where the 
technique yields unphysical results. In this limit the failure of 
the technique is due to the method neglecting analytic terms in 
the free energy, which in this case provide non-negligible con- 
tributions. In the spinless case, the variational results are un- 
physical when > 1/2 and 7 < 7c. As Fig.|5]illustrates, in 
this regime a exhibits two peaks per period. When gL > 1 this 
structure arises because the free energy is minimized when 
A = in the region of Vg between those peaks. As 7 de- 
creases, the width of this region increases and the peaks in 
a spread farther apart until A = everywhere and the peaks 
disappear entirely. This can be understood by noting that in 
the limit 7^1 the system essentially reduces to the single- 
impurity problem discussed in Section III for which A = is 
the only solution when > 1. In this limit the variational 
method thus predicts that Q = Vg- While this is consistent 
with renormalization group arguments when > 1, in the 
case of Fermi liquid leads (g^ = 1) one would expect to see 
small Fabry-Perot oscillations in the charge as the gate voltage 
varies. These Fabry-Perot oscillations have been observed in 
the conductance of single-walled nanotubes with near-perfect 
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Ohmic contacts-', and were also captured by a perturbative cal- 
culation of the charge in Section IVA. 

For gL < 1, the free energy is minimized by a nontrivial A 
for the entire range of gate voltages. This solution can be ob- 
tained analytically by first noting that A <C 1 when 7 ^ 1 so 
that the dominant contribution to the integral in /(A) comes 
from the region of small x where x < g. To approximate 
/(A), we therefore cut off the integral ?A x — g and expand 
K.±{2TTeLx) to first order in x. For Vg away from a half- 
integer, the square root term in Eq. ( I10> can be dropped, al- 
lowing determination of A. The resulting charge when Vq is 
away from a half-integer is 

Q - Vg oc rj^^-OL \ sin(7rVG)|| cos(7rVG)| • (Bl) 

When 1/2 < < 1, Eq. ( IBH predicts an unphysical double- 
peak structure in a similar to that shown in Figure|5jb), while 
for .gi < 1 /2, Eq. (IB It predicts oscillations in a with a single 
peak per period. An analysis of the charge at half-integer val- 
ues of Vg, however, reveals that the charge is discontinuous 
at these values when < 1/2. If we correct for this artifact, 
then the variational technique reproduces the expected Fabry- 
Perot oscillations in a when g^ < 1/2. 

In the spinful case, the variational technique yields unphys- 
ical results for the charge when 7' < 7^. Here, 7^ corresponds 



to the value of 7' below which the charge is given for all Vg 
by Eq. ( I18t . which has a nearly identical structure to Eq. Mil 
governing the charge in the spinless case. The behavior of 
(7 in this regime is therefore very similar to that in the spin- 
less case, and the method once again predicts an unphysical 
double-peak structure in a as depicted in Figure |9l The main 
difference here is that this failure occurs for all g^ since the ef- 
fective interaction parameter in the leads is [g^ + 1)/2 > 1/2. 
For (7l > 1, the region between the peaks in a where the free 
energy is minimized by A = increases as 7' is lowered. 
Eventually, A = everywhere and the peaks disappear alto- 
gether. 

In the limit 7' <C 1, the charge can be obtained analyti- 
cally by invoking similar approximations that were made in 
the spinless case. We again find only trivial solutions for A 
when (7l > 1 so that Q — 2Vg- For g^ < 1, however, the free 
energy is minimized by a nontrivial value of A, and the charge 
away from half-integer values of Vg is 

Q~2Vg'x 2r;7'T^|siii(7rVb)||cos(7rVG)|^- (B2) 

It follows from Eq. (IB2t that a exhibits an unphysical double- 
peak structure for all g^ < 1. When Vg equals a half-integer 
we find that A = is the only solution for any g^ < 1 so that 
the charge is a continuous function of Vg- 



